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ABSTRACT 

We present luminosity distributions for the X-ray sources detected with XMM-Newton in the bulge and 
disk of the Andromeda Galaxy (M31). The disk is clearly lacking the brighter sources which dominate 
X-ray emission from the bulge. This is the first convincing evidence for a difference between bulge 
and disk X-ray populations in M31. Our results are in good qualitative agreement with the luminosity 
distributions for low- and high-mass X-ray binaries recently obtained by Grimm et al.(2001) for our 
Galaxy. This confirms that X-ray population of the disk of M31 is dominated by fainter HMXB sources, 
while the bulge is populated with brighter LMXBs. 
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1. INTRODUCTION 

Multiple observations of external galaxies starting with 
Einstein have shown convincingly that the X-ray emission 
of normal spirals similar to our own M ilky Way is domi- 
nated by X-ray binaries ( Fabbiano 1995| ) . An X-ray binary 
contains either a neutron star or a black hole; it can be 
classified as low-mass (LMXB) or high-mass (HMXB) ac- 
cording to the spectral type of the companion star. It 
is well known that in our Galaxy most HMXBs lie in 
the Galactic plane, while the LMXBs have more spher- 
ical distribution with conspicuous concentration towards 
the Galactic Center. LMXBs and HMXBs demonstrate 
very similar X-ray spectra in some of their states, there- 
fore difficult to discriminate them in a single spectral ob- 
servation. Recently Grimm, Gilfanov and Sunyaev (2001) 
demonstrated the difference between X-ray luminosity dis- 
tributions of LMXBs and HMXBs in the Milky Way, pro- 
viding a method for discriminating the two populations in 
other galaxies. 

M31 (the Andromeda Galaxy) is of special importance 
for extragalactic astronomy because, at a distance of 760 
kpc (van der Bergh 2000), it is the closest spiral galaxy to 
our own. M31 is believed in many respects to be similar 
to the Milky Way. By observing M31 we can sample hun- 
dreds of X-ray sources at a nearly uniform distance and 
less obscured by interstellar gas and dust than those in 
the Galaxy due to the favorable orientation of the M31. 
It is therefore quite interesting to look at the luminosity 
distributions for various populations of X-ray sources in 
M31 and compare them with our Galaxy. However, earlier 
attempts to find the difference between luminosity distri- 
butions of X-r ay sources in the bulge and d isk of M31 were 
unsuccessful ( Lrinchicri fc Fabbiano 1991 ). 

In this Letter we report on the difference between the 
disk and bulge X-ray source populations of M31, as ob- 
served with XMM-Newton in the course of the most sen- 



sitive X-ray survey of a neighboring galaxy to date. 



2. OBSERVATIONS AND DATA ANALYSIS 

In the following analysis we use data from three XMM- 
Newton observations of the bulge and galactic disk of M31 
(see Table E], Fig. [IJa). The first observation of the cen- 
tral part of M 31 (M31 Core; r < 1 5') was performed on 
June 25, 2000 [Shirey et al. 200 1| , psborne et al . 2001 



Two regions of the northern part of the galactic disk 
(r > 15') were observed w ith XMM-Newton on January 5 
Trudolyubov et al. 200"2a] (M31 Northl Field) and January 



26, 2002 (M31 North2 Field) as a part of the Guaranteed 
Time Program (PI: K.O. Mason) (Fig. [[]). We use data 
from three European Photon Imaging Camera (EPIC) in- 
struments: two EPIC MOS detectors (|Turncr et al. 2001| ) 
and the EPIC PN detector (Ktrueder et al. 20011). In all 



observations the EPIC instruments were operated in the 
full window mode (30' diameter FOV) with the medium 
optical blocking filter. 

We reduced EPIC data with the XMM-Newton Science 
Analysis System (SAS v 5.2) 2 . Images in celestial coordi- 
nates with a pixel size of 2" have been accumulated in the 
0.3 - 12 keV energy band from both MOS and PN detec- 
tors. The resulting images for the EPIC-MOS1 instrument 
are shown in Figure 0&,c ; rf. Sources were detected with a 
program based on the wavelet decomposition algorithm, 
set at 4a detection threshold. We expect errors in the 
source position determination to be dominated by resid- 
ual systematic errors of order 2 — 5". 

We detected 230 point X-ray sources above a 4<7 level 
of significance. We cross-identified the detected X-ray 
objects with catalogs of globular cluster candidates and 
supernova remnants (Battistini et al. 1987, |Barmby &| 
Huchra 2001] , |Blair et al. 1981] , |Magnier et al. 1995| ), based 
on a positional coincidence of < 6" (Fig. p]6). Four sources 



1 Als o: Space Research Institute, Ru ssian Academy of Sciences, Profsoyuznaya 84/32, Moscow 117810, Russia 
2 See tittp:/ /xmm. vilspa.esa.es/usei 
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were identified with foreground objects and excluded from 
the analysis. Due to source confusion in the central part of 
the galaxy sources located within the central 1' were also 
excluded. 

To estimate energy fluxes and spectra, we used an ex- 
traction radius of ~ 20 — 60" (depending on the distance 
of the source from the telescope axis) and subtracted as 
background the spectrum of adjacent source-free regions. 
We corrected the count rates for the vignetting of the 
XMM telescope, based on the Current Calibration Files 
provided with SAS. The EPIC count rates were converted 
into energy fluxes using analytical fits to the spectra of 
brighter sources. Detailed spectral analysis for individual 
sources will be presented elsewhere. For faint sources we 
estimated energy fluxes with Web-PIMMS 3 assuming an 
absorbed simple power law spectral shape with photon in- 
dex a = 1.5 and an equivalent absorbing column of 7x 10 20 
cm~ 2 (Galactic value for the direction to M31). 

3. RESULTS AND DISCUSSION 

Altogether we detected 230 X-ray sources in the three 
fields of M31: 118 sources in the Core and 112 in the 
Northl and North2 fields. These include unidentified back- 
ground and foreground sources, as discussed below. The 
propert ies of the individua l sources are presented else- 



where (IShirey et al. 200 1[ |Osborne et al. 200l| , |fru 



dolyubov et al. 2002b). Here we discuss the source pop- 
ulations and their luminosity distributions. We have no- 
ticed that the brightest X-ray sources are concentrated 
towards the renter nf M31 , while the fainter ones are dis- 



trihntlcH much more uniformly (Fig 1 r,d\ Furthermore 

we have mentioned that all remaining bright sources in 
Northl and North2 fields can be identified with globular 
clusters (Fig. 16). Our data prove, for the first time, that 
M31 resembles our Galaxy in the distribution of its X-ray 
source populations ( i.e. LMXBs and HMXB) (see Fig.l in 
flGrimm et al. 200l|)). 



3.1. Luminosity distributions of X-ray sources 

We have built differential and cumulative luminosity dis- 
tributions of the detected X-ray sources for the central 
part of M31, and the combined Northl and North2 fields 
assuming a distance of 760 kpc (Fig. ||). In order to min- 
imize the contribution from non-disk X-ray population, 
the sources identified with globular clusters were excluded 
from the Northl/North2 sample. Taking into account the 
sensitivity as a function of off-axis distance, the different 
exposure times, and the effect of diffuse X-ray emission (in 
the M31 Core field), we estimate a flux completeness limit 
of our sample as ~ 10 36 ergs s _1 (indicated by the dotted 
line in Fig. 

Some of the objects in the M31 fields must be back- 
ground and foreground sources lying outside M31 (i.e. 
background AGN, and Galactic K and M stars ). Based 
on Chandra and XMM-Newton deep field results piacconj 



et al. 2001, Hasinger et al. 2001, one might expect to 



detect < 10 background AGN in each field with an appar 



ent luminosity around 10 ergs s . The total number of 
unidentifi ed foreground Galac tic stars could be ~ 5 sources 



per field Shirey et al. 2001. Correction for such back- 

3 See heasarc.gsfc.nasa.gov/Tools/w3pimms.html 



ground/foreground objects would flatten the faint end of 
the cumulative luminosity distribution (below ~ 3 x 10 36 
ergs s" 1 ), but would not significantly modify the higher 
luminosity part of the distribution. 

As it is clearly seen from both panels of Fig. ||, the pop- 
ulation of the northern disk regions of M31 is dominated 
by relatively faint objects with luminosities below ~ 10 37 
ergs s _1 . The overall shape of the cumulative luminos- 
ity distribution (Fig. ^6) can be presented with a simple 
power law function with an index of — 1.3 ± 0.2 4 for lumi- 
nosities above ~ 10 36 ergs s _1 . There is some indication 
for a flattening of the distribution at lower luminosities 
(Fig. |]&), which can be attributed to the incompleteness 
of our sample at the faint end. 

In contrast to the northern disk regions, bright X-ray 
sources contribute significantly to the population in the 
central region of M31 (Fig. ga). The shape of both cu- 
mulative and differential luminosity distributions for the 
central region of M31 {Blue histograms in Fig. 0) indicate 
the presence of a brake at around (1 — 2) x 10 ergs s _1 . 
We fit the unbinned cumulative distribution with a bro- 
ken power law (Fig. §6). For source luminosities between 
10 36 and 1.5 x 10 37 ergs s _1 , we obtain a power law index 
of —0.5 ± 0.1, while for the higher luminosities the power 
law index is — l^lg 2 . There is again some flattening of 
the cumulative distribution towards the faint source end 
(Fig. H&), probably caused by the incompleteness of our 
sample. These results for the M31 Core field are in gen- 
eral a greement with prev i ous studies of the central part of 
M31 (j Primini et al. 1993| , [Supper et al. 200l| , [Shirey et af 



200 1[ |Koiig et al. 2002| ) 



3.2. Comparison with previous results 

Studies of the luminosity distributions began as soon as 
individual X-ray sources in nearby galaxies were resolved 
(see Fabbiano 1995 for a review of early results). More 
sensitive surveys with ROSAT and Chandra significantly 
increased the number of sources detected (see Read, Pon- 
man & Strickland (1997) for ROSAT and Fabbiano (2001) 
for Chandra review). Most previous studies, however, dif- 
fer from our present work in two important aspects: (i) 
for most external galaxies, sensitivity limits were above or 
around 10 37 erg/s, while in our case we are able to study 
in detail population of sources between 5xl0 35 and 10 37 
erg/s; (ii) previous studies of luminosity distributio ns in 
M31 ( [rrinchieri k Fabbiano 1991| , iKong et al. 2002| ) were 
restricted to the central 17' — 20'. Since the inner region of 
M31 (r < 15') is dominated by the bulge and spheroid pop- 
ulations, the disk population was not well sampled in pre- 
vious surveys. Our results, however, can be directly com- 
pared to the luminosity distributions for Galactic sources 
obtained by Grimm et al. (2001). The comparison shows 
that the luminosity function which we obtain for the core 
region is close to the luminosity function of LMXBs in our 
Galaxy: both distributions have a characteristic "knee" 
at a fewxlO 37 ergs s _1 and can be described by a broken 
power law or a cut-off power law models. The luminos- 
ity distribution for the two northern fields of M31 resem- 
bles the luminosity function for Galactic HMXBs: both 
have a power law-like form with a slope steeper than that 



We fitted unbinned cumulative luminosity distributions using a maximum likelihood method. 
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of the corresponding M31 bulge/Galactic LMXB luminos- 
ity distributions; the M31 disk/Galactic HMXB distribu- 
tions terminate at much lower luminosities (<~ 10 37 ergs 
s _1 ) than M31 bulge/Galactic LMXB luminosity distri- 
butions (> 10 38 ergs s _1 ). We conclude therefore that 
the population of fainter X-ray sources we found in disk 
of M31 is probably dominated by the HMXBs, while the 
bulge region inside 15' is dominated by LMXBs by anal- 
ogy with the Milky Way. Our results confirm that, as in 
Milky Way and other nearby galaxies, the younger popula- 
tion of HMXBs in M31 is typically less luminous than the 
older populations of LMXBs and globular clusters (see also 
Hclfand & Moran 2001 for a discussion of HMXB luminosi- 
ties). In our study we see no evidence that the brightest 
sources in M31 are associated with young populations; in 
fact, we see the opposite effect. This is in contrast to So- 
ria & Wu (2002) , who invoke a population of bright young 
binaries to explain a large number of high luminosity ob- 
jects observed in the starburst galaxies. That effect can be 
probably attributed to the differences in the galaxy type 
and ranges of X-ray luminosities under study. 

4. CONCLUSIONS 

We analyzed XMM-Newton data for the three M31 fields 
and detected a total of 230 point-like sources. We built lu- 



minosity distributions separately for the bulge (inner 15') 
and disk regions. A striking difference between the bulge 
and disk populations has been demonstrated for the first 
time. We report on the discovery of lower luminosity pop- 
ulation in the disk of M31. Comparison of our results 
with the luminosity function obtained by Grimm et al. 
(2001) for the Milky Way allows us to conclude that lower- 
luminosity population in the disk is probably dominated 
by HMXBs, and that the luminosity distributions of X- 
ray sources in the bulge and disk of M31 are similar to the 
corresponding parts of our Galaxy. 
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Table 1 

XMM-Newton observations of M31 used in this analysis. 



Date, UT 


Tstart, UT 


Field 


Obs. ID 


RA (J2000) a 


Dec (J2000) a 


Exp.(MOS) 6 


Exp.(PN) 6 




(h:m:s) 






(h:m:s) 


(d:m:s) 


00 


(») 


25/06/2000 


10:44:42 


M31 Core 


0112570401 


00:42:43.0 


41:15:46.1 


34835 


31021 


05/01/2002 


06:28:31 


M31 Northl 


0109270701 


00:44:01.0 


41:35:57.0 


56810 


49568 


26/01/2002 


16:51:03 


M31 North2 


0109270301 


00:45:20.0 


41:56:09.0 


32193 


27982 



a - coordinates of the center of the field of view 



b - instrument exposure used in the analysis 
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Fig. 1. — Combined XMM-Newton/EPlC-MOSl images of the M31 Core, Northl and North2 Fields. Panel a is an optical image of M31 
from the Digitized Sky Survey (DSS) with EPIC-MOS1 field of view shown by solid black circles for each of the observations. The X-ray 
emission intensity is shown by logarithmic grey scale (black is the maximum) in Panels b,c,d representing enlar ged centra] (M31 C.nrp - r <r 1FV 

and northern pgir-tc Q f the galaxy (Northl and North2 fields). Positions of identified globular cluster candidates ( BalijatLaJ et " 1 1987 

2 ira 2001) are shown with small circles in Panel b. Identified supernova remnants (Blair et al 



Hue 



Barmby 



1981| , pMagmcr et al. 1995| ) are marked 
with small boxes in the same panel. In Panel c we circled all detected sources with luminosities greater than 10°' erg/s. Fainter sources are 
shown in Panel d. Some of the sources shown lie outside the EPIC-MOS1 sensitive area, but were detected with the EPIC-MOS2 or EPIC-PN 
detectors. It is evident that the brightest sources are concentrated towards the center of the galaxy, while fainter sources are distributed 
throughout the disk. A few bright sources in the disk regions (Northl and North2 fields) are identified with globular clusters, and hence 
belong to the non-disk subsystem of the galaxy. 
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Luminosity, xlO 35 (ergs/s) Luminosity, xlO 35 (ergs/s) 

Fig. 2. — Panel a: Differential luminosity distributions of the bulge (M31 Core) and disk (Northl and North2 regions combined) of M31. 
Panel b: Cumulative luminosity distributions for the same regions. A broken power law approximation to the luminosity distribution of the 
central part of the galaxy is shown with thin dashed line. A power law approximation to the luminosity distribution of the disk regions 
is shown with thin solid line. The indices of the power law fits to the luminosity distributions are also indicated. The vertical dotted line 
corresponds to the approximate completeness limit of our sample. 



